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Abstract 

We review recent studies demonstrating a nonuniversal (contin- 
uously variable) survival exponent for history-dependent random 
walks, and analyze a new example, the hard movable partial re- 
flector. These processes serve as a simplified models of infection in 
a medium with a history-dependent susceptibility, and for spread- 
ing in systems with an infinite number of absorbing configura- 
tions. The memory may take the form of a history-dependent 
step length, or be the result of a partial reflector whose position 
marks the maximum distance the walker has ventured from the 
origin. In each case, a process with memory is rendered Markovian 
by a suitable expansion of the state space. Asymptotic analysis 
of the probability generating function shows that, for large t, the 
survival probability decays as S{t) ~ where 5 varies with the 
parameters of the model. We report new results for a hard partial 
reflector, i.e., one that moves forward only when the walker does. 
When the walker tries to jump to the site R occupied by the re- 
flector, it is reflected back with probability r, and stays at R with 
probability 1 — r; only in the latter case does the reflector move 
(R R-l-1). For this model, 5 = 1/2(1 — r), and becomes arbi- 
trarily large as r approaches 1. This prediction is confirmed via 
iteration of the transfer matrix, which also reveals slowly-decaying 
corrections to scaling. 
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I. INTRODUCTION 



Random walks with absorbing or reflecting boundaries, or with memory, 
serve as important models in statistical physics, often admitting an exact anal- 
ysis. Among the many examples are equilibrium models for polymer adsorp- 
tion [1-4] and absorbing-state phase transitions [5]. Another motivation for 
the study of such problems is provided by the spreading of an epidemic in a 
medium with a long memory [6] . 

In addition to the intrinsic interest of such an infection with memory, our 
study is motivated by the spread of activity in models exhibiting an infinite 
number of absorbing configurations (IN AC), typified by the pair contact pro- 
cess [7,8]. Anomalies in critical spreading for IN AC, such as continuously vari- 
able critical exponents, have been traced to a long memory in the dynamics 
of the order parameter, p, due to coupling to an auxiliary field that remains 
frozen in regions where p — [8,9]. INAC appears to be particularly rele- 
vant to the transition to spatiotemporal chaos, as shown in a recent study of 
a coupled-map lattice with "laminar" and "turbulent" states, which revealed 
continuously variable spreading exponents [10]. 

Grassberger, Chate and Rousseau [6] proposed that spreading in INAC 
could be understood by studying a model with a unique absorbing configura- 
tion, hut in which the spreading rate of activity into previously inactive regions 
is different than for revisiting a region that has already been active. In light 
of the anomalies found in spreading in models with INAC or with a memory, 
we are interested in studying the effect of such a memory on the scaling be- 
havior in a model whose asymptotic behavior can be determined exactly. Of 
particular interest is the survival probability S{t) (i.e., not to have fallen into 
the absorbing state up to time t). 

In the present work, we review previous results on survival of random walks 
with memory, and analyze the asymptotic behavior of a random walk subject to 
a hard movable refiector. On the basis of an exact solution for the probability 
generating function, we obtain the decay exponent S. 

The balance of this paper is organized as follows. Section II reviews previ- 
ous results on variable survival exponents for random walks with memory. In 
Sec. Ill we define the hard reflector model and obtain a formal solution for 
the generating function. An asymptotic analysis of this function is presented 
in Sec. IV, leading to an expression for the decay exponent S in terms of the 
reflection probability r. In Sec. V we present exact numerical results for fl- 
nite times (from iteration of the probability transfer matrix) that complement 
and extend the asymptotic analysis. Sec. VI contains a brief summary and 
discussion. 
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II. VARIABLE SURVIVAL EXPONENTS IN RANDOM WALKS 

WITH MEMORY 



Relatively simple random walk problems often serve as reduced examples 
of much more complicated many-body phemomena. So it is with phase tran- 
sitions between an active and an absorbing state. In these systems [5,11] the 
stationary state of a Markov process exhibits (in the infinite-size limit) a phase 
transition from a frozen, inactive regime to one with sustained activity, as a 
control parameter is varied. Broadly speaking, the control parameter repre- 
sents the reproduction rate of activity (A — > 2A) relative to its extinction 
(A — > 0). The simplest examples are the contact process [12] and directed 
percolation (DP). 

Consider now, instead of the stationary state, the spread of activity from a 
localized source, in an infinite system. In the subcritical regime (for which the 
only stationary state is the inactive, absorbing one) , the initial activity decays 
(typically, exponentially fast), while in the supercritical regime there is a finite 
probability for it to spread indefinitely. Just at the critical point, one finds a 
scale- invariant evolution: the survival probability S{t), the integrated activity 
n{t), and the mean-square distance, R'^{t), of the activity from the initial seed, 
all follow asymptotic power laws. 

The question of survival arises naturally in the context of a random walk 
in the presence of an absorbing boundary. The survival probability S{t) is 
the probability never to have visited an absorbing boundary until time t. The 
simplest example is a one-dimensional random walk Xt (in discrete time) on 
the non- negative integers, with the origin absorbing, and Xq > 1. Let the 
walker jump, at each time step, to the right {xt+i = Xt + 1) with probability 
p, and the left {xt+i = Xt — 1) with probability q = 1 — p. U p < q then the 
survival probability decays exponentially, while for p > q it approaches (again 
exponentially) a nonzero value, so that the walker has a finite probability to 
escape to infinity. In the absence of drift {p — q — 1/2), S{t) ~ for large 
t; associated with this power-law decay is an infinite mean lifetime. In the 
analogy with an absorbing-state phase transition, p = 1/2 evidently marks the 
transition, with extinction certain for p < 1/2, and a finite asymptotic survival 
probability for p > 1/2. The same qualitative situation holds in the contact 
process, starting for example from a single active site [5] . 

The analogy is in fact exact for the rather special case of compact di- 
rected percolation (CDP), in which active regions are delimited by independent 
random walks that annihilate on contact. (CDP is a particular limit of the 
Domany-Kinzel model [13].) In this case 'drift' corresponds to a tendency of 
the walkers at the boundaries of an active region to approach, or separate from, 
one another; the critical point corresponds to zero drift, or unbiased random 
walkers {p = 1/2). Initializing CDP with a finite interval (say, 1, 2,...,m) of 
active sites, is equivalent to placing random walkers {xt and yt) at and m+1. 
The active region at any subsequent time corresponds to the interval between 
the walkers Xf and y^; the process ends when the two meet. To make the anal- 
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ogy between CDP and a random walk with the origin absorbing complete, we 
may fix Xf = for all times, so that only the right frontier of the active region 
is free to fluctuate, while the left frontier is pinned, as it were, at a wall. 

Given the connection with phase transitions, we shall think of the power law 
for the survival probabilty as defining a critical exponent, and write S{t) ~ t~^. 
The main interest of the examples discussed in this paper is that the exponent 
S can be shown to vary continuously as a function of a parameter. This in turn 
may help to understand scaling in more complex examples, such as the pair 
contact process [7] , for which a variable survival exponent has been reported 
numerically, but cannot be established rigorously. 

Random walk models exhibiting variable survival exponents fall in two 
classes. In one, the position of the absorbing boundary is a given (determin- 
istic) function of time. A random walk in the presence of such a boundary 
defines a nonstationary stochastic process. The second class (which is our 
principal interest here) involves memory, either in the form of a refiector that 
moves when it encounters the walker, or of a history-dependent step length. 

We begin with a brief review of the first class. In a highly readable paper, 
Krapivsky and Redner [14] considered what happens when a random walker on 
the line is subject to two absorbing boundaries, R{t) and L{t), which are pre- 
scribed functions of time (initially the walker is at the origin) . The absorbers 
are initially near the origin, and follow R{t) = a + {At)", L(t) = —R{t). For 
a > 1/2, the absorbers rapidly leave the region explored by the walker, which 
therefore enjoys a finite probability of survival as t ^ oo. If the absorbers are 
stationary we of course expect S{t) to decay exponentially; but if their motion 
is characterized by0<ci;<l/2, this changes to a stretched exponential decay, 
S ~ exp[— i-^"^"]. For a = 1/2, the 'safe' region expands with the same power 
law as the region explored by the walker. The result is power-law decay of the 
survival probability, with a variable exponent which depends on A. 

The problem studied by Krapivsky and Redner may be pictured as a ran- 
dom walk confined to a parabola in the t — x plane, whose equation (for t > 0) 
is X — ±{At)°'. When a = 1/2, as noted, there is power-law decay of S(t), with 
a nonuniversal exponent. Similar conclusions apply for DP and for directed 
self- avoiding walks [15], and for CDP [16], when these processes are confined 
to a fixed parabola. 

We now turn to studies of a random walk subject to some special condition 
when it enters virgin territory, i.e., when it attempts to visit a site for the first 
time. Here there is no fixed confining boundary; the condition depends on the 
history of the walk. But since the region explored by the walker grows oc t^/^, 
it effectively creates its own parabola. 

It was recently shown [17] that an unbiased random walk on the non- 
negative integers, with the origin absorbing, exhibits a continuously variable 
exponent S when subject to a mobile, partial refiector. The latter is initially 
one site to the right of the walker. Each time the walker steps onto the site 
occupied by the refiector, it is refiected one step to the left with probability 
r (it remains at its new location with probability 1 — r) ; in either case, the 
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reflector is pushed forward one site in this encounter. The survival exponent 
S = (l + r)/2 in this process [17]. Since the reflector effectively records the 
span of the walk (i.e., the rightmost site yet visited), its interaction with the 
walker represents a memory. We shall refer to this model as the soft reflector, 
to highlight the fact that the refector obhgingly moves forward, even if the 
walker is reflected back. Note that in this case S never exceeds unity. 

The analysis of the random walk with soft reflector was subsequently ex- 
tended to compact directed percolation [18]. The active region initially consists 
of a single site (the origin), and, as already noted, is bounded by a pair of in- 
dependent, unbiased random walkers, originally ai x — and x — 1. The two 
walkers are subject to movable partial ("soft") reflectors, such that the walker 
on the right is reflected toward the left and vice-versa. The results for the 
survival exponent are qualitatively similar to those for the single walker, but 5 
now varies between 1/2 and 1.160 as the reflection probability r varies between 
zero and one. The results (coming in this case from iteration of the transition 
matrix, rather than from an asymptotic analysis of the generating function), 
are well flt by the simple expression 6 = l/2 + 2r/3; small but signiflcant 
deviations from this simple formula are found, however. CDP with reflectors 
has so far defled exact analysis, and the reason for the speciflc value S = 1.160 
for r = 1 is not understood. 

Most recently, the methods developed in Ref. [17] were applied to a one- 
dimensional random walk with memory of a different form: if the target site 
X lies in the region that has been visited before (that is, if x itself has been 
visited, or lies between two sites that have been visited), then the step length 
is v; otherwise the step length is n. In this case one finds 5 — v/2n [19]. Thus 
S can take any rational value between zero and infinity. 

With this background we may now describe the problem to be analyzed 
here as a random walk subject to a hard partial, movable refiector. This is 
because the refiector now moves forward if and only if the walker succeeds in 
occupying the new position; when the walker is reflected, the reflector main- 
tains its position. This, as will be shown, can lead to much larger values of 5 
than in the soft reflector case. 

Before embarking on the technical discussion, we summarize our approach, 
as developed in Refs. [17] and [19]. After formulating the problem, we enlarge 
the state space so that the process becomes Markovian in the expanded repre- 
sentation [20] . We then write down the equation of motion for the probability 
distribution and its associated boundary and initial conditions. Since these 
are discrete models, the equation of motion corresponds to a set of difference 
equations, flrst-order in time, and second-order in space. It is convenient to 
eliminate the time variable by passing to a generating function P{z) (effec- 
tively, a discrete Laplace transform). Using separation of variables, we obtain 
a formal solution for the generating function. Finally, the asymptotic long- 
time behavior is found by studying the generating function for the survival 
probability in the limit z ^ 1. 
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III. MODEL 



Consider an unbiased, discrete-time random walk on the nonnegative inte- 
gers, with the origin absorbing. We denote the position of the walker at time t 
by Xt, with xq — 1. The movement of the walker is affected by the presence of 
a movable partial reflector, whose position is denoted by Rt, with Rq = 2. At 
each time step the walker hops from its current position Xt to cither Xt + 1 or 
Xt — 1 with probabilities of 1/2. If, however, xt + l = Rt, the walker is reflected 
back to Xt with probability r, and remains at a:;^ -|- 1 with probability r = 1— r; 
in the latter case the reflector simultaneously moves to Rt + 1. Summarizing, 
the transition probabilities for the walker are 

_{xt-l , w.p. 1/2 
""'^'-{xt + l ,w.p. 1/2 

in case Xt < Rt — 2. When Xt = Rt — 1, we have instead 

{Xt-1 , w.p. 1/2 
Xt + l , w.p. r/2 (2) 
Xt , w.p. r/2 

The position of the reflector at any moment is given by -Rt = 1 + maxt'Ktixf}. 

Although the process Xt is non-Markovian (since the transition probability 
into a given site depends on whether it has been visited previously), we can 
define a Markov process by expanding the state space to include the variable 
Ut = Rt — 1 — meiXt'<t{xt'}- The state space E al? is given by by 

E = {{x,y)eI?:x>Q, y > 1, x<y}, 

as represented in Fig. 1. 

The probability distribution P{x,y,t) follows the evolution equation 

P{x, y,t+l)^ ^P{x + 1, y, t) + ^P{x -l,y,t), ior x < y - 1, (3) 

with P{x,y,0) = Eq- (3) is subject to two boundary conditions. The 

first is the absorbing condition for a; < 

P{x,y,t)^0, forx<0. (4) 

The second applies along the diagonal x — y. Defining D{y,t) = P{y,y,t), we 
have 

D{y, t+l) = \p{y-l, ^) + ^ D{y-l, t) + t), for y>2. (5) 

For y —1 the equation is simply t-\-l) — (r/2)£)(l, t), and since D(l, 0) = 
1, one has D{1, t) — (r/2)*. Finally, ior x — y — 1, 

P(y- 1, y,t+l) = \p{.y-2, y, t) + ]^D{y, t), for a; < y - 1, (6) 
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We next introduce the generating function: 



P(x,y,^) = ^P(a;,y,i)z*. (7) 
t=o 

Multiplying Eqs. (3), (6) by 2;*, summing over t and shifting the sum index 
where necessary, one finds that the generating function satisfies 

-P{x,y) = lp{x + l,y) + lp{x-l,y), iorx<y-2 (8) 
z z z 

-P{y-h y) = \D{y) + ^^(z/-2, y), for x = y - l, (9) 

z z z 

(we drop the argument z for brevity), where D[y) is defined by an expression 
analogous to Eq. (7). The initial condition implies D{1) = (1 — zr/2)~^; the 
boundary conditions are 

P(0,y) = 0, (10) 

and 

-b{y) = ]-P{y-l,y) + "^biy-l) + ^^'(y), for y > 2. (11) 

z z z z 

Eq. (8) relates P at different values of for the same y. Specifically, on 
the interior of each line of constant y, P satisfies a diffusion equation, with 
a source at x = y and a sink at x = 0. It is therefore natural to attempt 
separation of variables, 

P{x,y)=A{x)B{y). (12) 
Inserting this expression in Eq. (8) one obtains 

-A{x) - \a{x - 1) - \a{x + 1) = , (13) 
z z z 

with A{Qi) = 0. The solution satisfying this boundary condition is 

A{x) = - A-^ , (14) 

with 



z \ z"^ 



Our next task is to determine P(y); for this we require a relation between 
generating functions with different arguments y. Relations of this kind arise 
along the diagonal, but involve the function -D(y), which we proceed to elimi- 
nate. Combining Eqs. (9) and (11), we find 
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Q{z)D{y) = -D{y-1) + jP{y-2,y) 



(16) 



where 



Equation (11) may also be written as 

(l-^) D{y)-'-LD^y-l)^'-P{y-l,y), 



(17) 



(18) 



We now multiply Eq. (16) for y — 1 by P = zr/{2 — zr) and subtract the result 
from the same equation for y. This yields 



Q 



D{y)-(3D{y-l) 



zr 

T 



b{y-l) - (3b{y-2) 



^From Eq. (18) we have 

b{y)-0{y-l) 



2- zr 

allowing us to eliminate 1) from Eq. (19): 



P{y-2,y)-pP{y-3,y-l) 
(19) 

(20) 



4Q 



z{2 - zr) 



P{y-l,y)-P{y-2,y)=(3 



'-P{y-2,y-l)-P{y-3,y-l) 



Inserting Eq. (12) one readily finds a recursion relation for B: 
B{y) _ zf[2A{y-2) - zA{y-?>)\ 



B{y-l) (4 - 2zr - z'^)A{y-l) - z{2 - zr)A{y-2) 



(21) 



(22) 



Given B{1) = i)(l)/i(l) with b{l) and A{y) as found above, Eqs. (12), 
(14), and (22) represent a complete formal solution for the generating function 
P{x,y). 



IV. ASYMPTOTIC ANALYSIS 

Our goal is to find the survival probability S(t) for large t. This can be 
found analysing the associated generating function. 



S{z) = J2S{t)^ 



(23) 



t=o 
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in the limit z ^ 1. Specifically, if S{t) ~ t~^, then the radius of convergence of 
S{z) is \z\ = 1, and the singular behavior of the generating function as z — > 1 
determines S. Indeed, in this case, with z = 1 — e, we have 

oo 

t=i 

/oo 
dt r'^^cxp[-t| ln(l - e)|] 

~ e'^-^r(l - S) , (24) 

so that the scaling exponent 5 can be read off from the power-law dependence 
ofjS'one^l — zase— >0. This simplifies considerably the determination of 
the long-time asymptotic behavior of S{t). 

The generating function S has two contributions, coming from the "inte- 
rior" {x < y) and the diagonal: 

oo 3/— 1 oo 
j/=l x=l y=l 

= Sp + Sd. (25) 
Using Eq. (20) one readily shows that 

oo 

'^D = :r^$^P(2/-l,2/)<-Sp, (26) 
^ ^ y=l 

so that it suffices to analyze the behavior of Sp. 
Consider 

oo y-l 

^P = 5]S(y)J]i(a;). (27) 

y=l x=l 

The sum over x can be evaluated as 

x=l 

^isinh^^ (28) 
A 2 ^ ^ 

where in the last step we used A = InA ~ \/2e as e ^ 0. 

To evaluate Sp we also require an expression for B{y) in the limit e ^ 0; 
this can be obtained from Eq. (22). We begin by setting all explicit factors of z 
equal to unity, since the 0(e) corrections thereby discarded do not contribute to 
the singular behavior of Sp. The singular contributions in fact originate from 
the functions A, through their dependence on A. Writing A{y) — 2sinhAy, we 
therefore have 
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B{y) ^ f[2sinhA(y-2)-smhA(y-3)] ^^g) 



B{y-l) (3-2r)sinhA(y-l)-(2-r)smhA(y-2) ' 

Using the identity sinh(a + h) — sinhacoshft + sinhftcosha, and neglecting 
terms C(A^), we obtain 

Biy) tanhAy-A , , 

vy; y (3Q) 



B{y-1) tanhAy + A/r ' 

For y < yo = [2/r] + 1 (here [...] denotes the integer part of its argument), 
we can write, for small A 

Biy)=mi{^^-CB{1), (31) 

k=2 ' 

where C depends on r and y but is independent of A. We shall in fact discard 
the contribution due to y < yo in Sp. The reason is that the contribution to 
the survival probability from any fixed, finite set of transient states must decay 
exponentially at long times, and so will not affect our result for the scaling 
exponent. 

Noting that ^(1) = £)(l)/i(l) oc 1/A, we have, for y > yo, 

, C A B(k) , , 

Biy) = -r M . ^ ^ , 32 
A 5(fc-l) ^ ' 

k=yo ^ I 

where C is a constant. Since all terms have k >2/f, we may use Eq. (30) to 
write, with 0^ = tanhA/c, 

BiVo) ty, l + A/(r0.) 

^-fi + =VE:^- (33) 

Approximating the sum by an integral we find 

+ (34) 

Now, inserting Eqs. (28) and (34) in Eq. (27), the generating function for 
e — > is: 

q Ai/f-iV^ sinh^Ay/2 

''^^ ^sinh^-V^A,' (35) 

where '~' denotes asymptotic proportionality as e — > 0. Approximating, as 
before, the sum by an integral, we have 
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Sp ~ A^/^-^ / / . (36) 

Since 1 + 1/f > 2, the integral converges at its upper limit. For r < 1/2, 
1 + 1/f < 3 and the integral remains finite as A — > 0. Then 

so that the survival probability decays as S{t) ~ with S — l/2f. For 
r — 1/2, the prefactor in Eq. (36) is independent of A and 



poo 



— ~-ln(l-^). (38) 



u 



Expanding the logarithm, we find S{z) ~ Xln-^"/'"'' yielding directly S{t) ~ 
t"^. Finally, when r > 1/2, 1 + 1/f > 3 and the integral in Eq. (36) contains 
two principal contributions: one finite (due to the interval from say, 1, to infin- 
ity), the other arising from the lower limit, and diverging as A^~^/^. Combined 
with the prefactor oc A^^^~^ however, the latter contribution is nonsingular, 
while the former is again proportional to A^/*""^. 

Summarizing, the asymptotic survival probability decays as a power law, 

s{t) ~ r^/^ , (39) 

which is the result we set out to prove. 



V. NUMERICAL RESULTS 

The foregoing analysis provides the t — > oo scaling behavior of the survival 
probability, but does not indicate the rate of convergence to the asymptotic 
power law. To determine how the corrections to scaling decay, we iterate the 
discrete time evolution equation for P{x,y,t) directly. In Fig. 2 we show the 
decay of S{t) for refiection probability r = 0.85, corresponding to 5 — 10/3. 
For very late times, the graph indeed approaches a power law with the expected 
exponent. The approach is, however, extremely slow. 

Since our asymptotic analysis only retains the leading dominant term in 
the long-time behavior of S{t), we have no specific information on correction 
to scaling terms. It is easy to see, nonetheless, that corrections oc t~^^'^ will 
be generated, since A ~ ^/2e + 0{e). In fact, we are able to fit the long-time 
evolution of the survival probability by adding a suitable 0{t~^^'^) term to the 
power law, but further terms (~ t^^), etc.) are required for intermediate times. 

We have found a particularly simple transformation of variable that appears 
to take the dominant correction to scaling into account. It consists in defining 
a shifted time variable 

T^t + bt^/^, (40) 
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with parameter b adjusted to make the graph of S versus T (on log scales) as 
linear as possible. Fig. 3 shows the survival probability data for r = 1/2, 2/3, 
0.75 and 0.85 versus T (the corresponding b vahics are 1.754, 4.167, 7.042, and 
16.67). In each case the numerical data (points) follow the modified power 
law, 

S{T) = AT-^, (41) 

to very high precision. [Here A is an amplitude determined by extrapolating 
T^S{T) to T ^ oo.] The numerical data appear to converge rapidly (faster 
than a power law) to the fit. While this 'shifted time' analysis is for the mo- 
ment without theoretical basis, it clearly confirms the asymptotic power laws 
found analytically, and suggests a simple form for describing slowly decaying 
corrections to scaling. 

The slowly decaying correction to scaling would likely frustrate efforts to 
extract the correct long-time behavior from simulations. Looking at Fig. 2, 
we see that the asymptotic power law is barely evident when S{t) has decayed 
to e"^^. To obtain even marginally useful simulation data in this situation 
we would need to perform > lOe^^ ~ 3 x 10^ independent reahzations of the 
process, extending to a maximum time of about 2000 steps. This is feasi- 
ble for a simple random walk, but becomes a computational challenge for a 
many-particle system. Thus, if lattice models such as the PCP behave in a 
manner analogous to what is found for the random walk with a hard reflector, 
it will be very difficult to confirm power-law scaling in simulations. Data for 
limited times (or limited samples) may well give the impression of faster than 
exponential decay of S{t). 

VI. DISCUSSION 

We have reviewed examples of confined random walks, and random walks 
with memory, that lead to a continuously variable scaling exponent for the 
survival probability, and analysed in detail the 'hard reflector' case. The latter 
problem appears to be particularly relevant to spreading in the pair contact 
process, since modification of the background density of isolated particles can 
only occur when activity invades a previously inactive region [7]. The strong 
correction to scaling found numerically for the hard reflector is reminiscent of 
the slow convergence (interpreted as faster than power-law decay in Ref. [6]), 
found in spreading studies of the PCP. 

Several interesting issues remain open. First, the nature of correction to 
scaling terms needs to be investigated using a more complete asymptotic ex- 
pansion of the generating function. Second, one would like to understand the 
exponent vahics for CDF (obtained numerically in Ref. [18]) on the basis of 
the generating function approach. Finally, extension of any of the models dis- 
cussed here to two or more dimensions promises to be a difficult but potentially 
fascinating challenge. 
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Figure Captions 

Fig. 1. Random walk subject to hard reflector: transitions in the x-y plane. 

Fig. 2. Decay of survival probability S{t) for reflection probability r = 0.85 
(solid curve); the slope of the straight line is -10/3. 

Fig. 3. Survival probabihty 5" as a function of the shifted time variable T, for 
reflection probabilities r = 1/2, 2/3. 3/4 and 0.85 (data points); the straight 
lines have slopes of -1, -3/2, -2, and -10/3. 



15 




1 2 3 4 5 6 7 




"6 
■12 
■18 
■24 
■30 



